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Abstract 


In this paper, a transient two-phase non-isothermal PEM fuel cell model has been developed based on the previously established two-phase 
mixed-domain approach. This model is capable of solving two-phase flow and heat transfer processes simultaneously and has been applied herein 
for two-dimensional time-accurate simulations to fully examine the effects of liquid water transport and heat transfer phenomena on the transient 
responses of a PEM fuel cell undergoing a step change of cell voltage, with and without condensation/evaporation interfaces. The present numerical 
results show that under isothermal two-phase conditions, the presence of liquid water in the porous materials increases the current density over-shoot 
and under-shoot, while under the non-isothermal two-phase conditions, the heat transfer process significantly increases the transient response time. 
The present studies also indicate that proper consideration of the liquid droplet coverage at the GDL/GC interface results in the increased liquid 
saturation values inside the porous materials and consequently the drastically increased over-shoot and under-shoot of the current density. In fact, 
the transient characteristics of the interfacial liquid droplet coverage could exert influences on not only the magnitude but also the time of the 


transient response process. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Numerical modeling and simulation of PEM fuel cells is cru- 
cial for revealing the underlying physics and facilitating cell 
design and optimization. Many PEM fuel cell models have been 
developed to achieve these goals, including the single-phase or 
pseudo single-phase PEM fuel cell models [1-13], which either 
do not account for liquid water transport phenomena or simplify 
its treatment, and the two-phase models [14—26], which handle 
the liquid water transport and the related phenomena explicitly. 
The majority of these models, however, consider only steady- 
state operation conditions of PEM fuel cells. Since a PEM fuel 
cell would experience start-up and shut-down processes and 
frequent load changes for automotive applications, a numeri- 
cal model capable of investigating transient responses is thus 
needed. 
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A transient PEM fuel cell model was developed in the work 
of Um et al. [1], but the main studies were not focused on the 
transient analysis, and thus only a preliminary result concerning 
the cell dynamic response with a step change of cell voltage was 
presented. This work was further extended by Wang and Wang 
[27] into a transient three-dimensional single-phase isothermal 
PEM fuel cell model. They estimated three time constants con- 
cerning electrochemical double-layer charging/discharging, gas 
transport, and membrane hydration/dehydration processes, and 
concluded that the double-layer charging/discharging process 
was very fast and could thus be safely neglected in transient 
analyses of PEM fuel cell operations. They also conducted 
extensive numerical simulations with step changes of cell volt- 
age and inlet humidity and discussed the dynamic physics of 
the transient phenomena. Wang and Wang [28] later studied the 
dynamic responses of a PEM fuel cell undergoing a step change 
of current based on the work of Meng and Wang [5,29] with a 
proper account of the electron transport. Shimpalee et al. [30,31] 
numerically simulated the transient responses of a PEM fuel cell 
subjected to variable load changes based on a three-dimensional 
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Nomenclature 


a 


Y'S 
£ 


EQNS ® 


BSS 


water activity or stoichiometry coefficient 
molar concentration (mol m~?) 
constant-pressure heat capacity (Jkg~! K7!) 
mass diffusivity (m s7!) 

water content diffusivity (mol m`! s7!) 
equivalent weight of the membrane (kg mol!) 
Faraday constant (96,487 C mol`!) 
transfer current density (A m~?) 

thermal conductivity (W m`! K7!) 
permeability (m7) 

electro-osmotic drag coefficient 
gas-phase pressure (Pa) 

capillary pressure (Pa) 

universal gas constant (J mol`! K7!) 
liquid saturation 

source term 

time (s) 

temperature (K) 

gas-phase velocity (m s7!) 

open-circuit potential (V) 

cell voltage (V) 

molecular weight (kg mol~!) 


Greek letters 


E porosity 

Em fraction of the membrane phase in the catalyst 
layer 

n over-potential (V) 

K proton conductivity (S m7!) 

À water content 

u viscosity (kg m7! s7!) 

Oc contact angle 

p gaseous density (kg m~?) 

o electronic conductivity (S m~ 1) or surface tension 
Nm’) 

T viscous stress tensor 

b phase potential (V) 

Superscripts 

cl catalyst layer 

eff effective value 

l liquid phase 

sat saturation value 

v vapor phase 

Subscripts 

cl catalyst layer 

e electrolyte or energy 

g gaseous phase 

i species 

l liquid 

m membrane 

s electron or solid phase 

sat saturation value 

w water 


single-phase isothermal PEM fuel cell model. The simulations 
were conducted for a PEM fuel cell with a serpentine flow-field 
and a 10cm? reactive area. They presented numerical results 
under excess, normal, and minimal fuel and air supplies. 
Recently, Wu et al. [32] developed a transient two-dimensional 
single-phase non-isothermal model of PEM fuel cells, and with 
an inclusion of the heat transfer equation, they were able to 
consider four transient processes, namely the electrochemical 
double-layer charging/discharging, species transport, mem- 
brane hydration/dehydration, and heat transfer processes. Based 
on their numerical results, they concluded that the heat transfer 
process could exert significant influences on fuel cell dynamic 
responses. 

The transient PEM fuel cell models briefly reviewed in 
the early section are all based on single-phase or pseudo 
single-phase simplifications. Wu et al. [33] have attempted to 
consider the liquid water effect by defining a fixed liquid sat- 
uration of 10% in the porous materials in a PEM fuel cell 
based on the two-phase results provided in the review paper 
of Wang [34]. Natarajan and Nguyen [21] developed a tran- 
sient two-dimensional two-phase model for the cathode of a 
PEM fuel cell, and they concluded that liquid water trans- 
port would prolong the cell response time, especially under 
the current-collecting land. This work only considered the 
cathode side under isothermal conditions. Since the inclusion 
of the temperature effect is crucial for correctly simulating 
the condensation/evaporation phenomena in two-phase flows, 
Shah et al. [35] recently developed a transient non-isothermal 
model for a PEM fuel cell. They presented numerical results 
in the form of potential sweeps and were able to quantita- 
tively predict the hysteresis phenomenon often observed in PEM 
fuel cell experiments. However, the model is one-dimensional 
and thus could not be applied to fully investigate the com- 
plex multi-dimensional physics in practical PEM fuel cell 
operations. Song et al. [36] also developed a transient one- 
dimensional non-isothermal two-phase PEM fuel cell model 
to investigate transient liquid water transport in the cathode 
GDL. The phase change phenomenon and parameters affecting 
transient and steady-state liquid water transport were dis- 
cussed. 

In order to fully investigate the dynamic responses of a 
PEM fuel cell under practical operation conditions and further 
enhance fundamental understandings of the intricate interac- 
tions of thermal and water managements during the cell transient 
operations, a transient two-phase non-isothermal PEM fuel cell 
model is developed in the present paper based on a previously 
established two-phase mixed-domain approach [26]. The model 
is capable of solving two-phase flow and heat transfer pro- 
cesses simultaneously with a proper consideration of the effect 
of liquid droplet coverage at the gas diffusion layer (GDL) 
and the gas channel (GC) interface. The model will be applied 
herein for two-dimensional time-accurate simulations in a cross- 
section perpendicular to the flow direction so that the effects 
of liquid water transport and heat transfer phenomena on the 
transient responses of a PEM fuel cell, with and without con- 
densation/evaporation interfaces, could be fully examined and 
clearly presented. 
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2. Theoretical formulation 


The complete conservation equations in their transient forms 
are developed in this paper based on a previously established 
steady-state multi-dimensional two-phase mixed-domain PEM 
fuel cell model [26]. First, the transient conservation equations 
of mass, momentum, species concentrations in the gaseous phase 
are presented. 


Mass conservation: 


d[ed — s 7 
GLI — OPN oy ety so (1) 
ot 
Momentum conservation: 
1 a(pu) 1 wes 
V- =-V V. S 
e(1 —s) ot e2(1—s) (puu) p+ T+ Su 
(2) 
Species conservation: 
aled — sci 2 
Be ey oa av aes 3) 


ot 


In Eq. (3), based on the mixed-domain approach [12,13], the 
water vapor concentration is solved only in the gas channels, gas 
diffusion layers, and catalyst layers on both the anode and cath- 
ode sides. In the two catalyst layers, the dissolved water phase 
(water in the membrane phase) is assumed to be in thermo- 
dynamic phase equilibrium with water vapor, and its transport 
process is combined into the water vapor transport equation 
using the following water diffusivity [12,13]: 
pol = l5 pele + el p, Rut & (4) 

Psat da 
In addition, the effective porosity, sf, in the water transport 
equation should be modified as [27,32] 


eL py Em RuT dà 


(5) 
1— sS Psa da 

In other species transport equations, the effective porosity 

remains as 


efe (6) 


Considering the liquid water effect, the effective gaseous 
species diffusion coefficients should be further modified as 


DË = Di — s)!> (7) 


In the present two-phase model, water produced in the 
cathode catalyst layer is assumed to be in vapor phase or 
dissolved phase [26], as they are at thermodynamic equilib- 
rium. Although the thermodynamic equilibrium condition is a 
common assumption also made in many prior transient numer- 
ical models [27,28,32,33], its validity in transient simulations 
requires further verification. 

In this transient two-phase model, liquid water transport is 
taken into account using the traditional two-fluid method. As 
discussed in [26], we will ensure the states of water vapor (or 


the dissolved water phase in the catalyst layers) and liquid water 
closely approach their thermodynamic equilibrium conditions 
by choosing appropriate condensation and evaporation rate coef- 
ficients. It should also be emphasized that the validity of this 
assumption needs further investigation for transient simulations. 
The transient conservation equations of liquid water transport 
are presented next. 
Liquid mass conservation: 
d(Epis) a 
“a + V- (ou) = SviWw (8) 
where an expression for the condensation/evaporation rate can 
be found in [26]. In the porous materials, based on the Darcy’s 
law, the liquid water velocity is derived as 
= KaK 1 


u= = Vp (9) 
H 


where the liquid pressure is defined as [14] 
P =p- Pe (10) 


Eqs. (8)—(10) can be combined to produce a transient conser- 
vation equation for the liquid saturation, which is 


ð KaK a KaK 
Cad y, [2 nl Pevs < E rl v| _ 


ot 


The transient water content conservation equation inside the 
membrane is in the following form: 


0 / Pmà 


- (2) =V: (DVA) +S a2) 


The transient conservation equation of energy is derived as 
3 : > 3 
aPC os TIF V (PC pT) + V - (Cp 1T) 


=V-(kK#VT) + Sr (13) 


where the effective parameter, (oC Pap can be determined in 


the following general form: 
(oC p)™ = e(1 — 8)pCp + €5(0Cp), + (1 — € — Em\(0Cp), 
+ Em(PCp) (14) 


The related parameters are provided in Table 1. The effective 
thermal conductivity of each cell layer is approximated as a con- 
stant in the present numerical studies and is also listed in Table 1. 
It should be noted that, in this paper, specific heat capacity, Cp, of 
the membrane material has been substituted with that of PTFE, 
since it is not available in the open literature, as also discussed 
in [32]. 

Finally, the conservation equations of proton and electron 
transport are derived as 


Proton transport: 


V- (°F Voe) + Se = 0 (15) 
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Table 1 
Physicochemical parameters 


Anode volumetric exchange current density, ajo (A m~) 1.0E+9 
Cathode volumetric exchange current density, ajo (A m~?) 1.0E+4 
Reference hydrogen concentration, Cy, (mol m7?) 40 
Reference oxygen concentration, Co, (mol m~?) 40 
Anode transfer coefficients a= =l 
Cathode transfer coefficient ae=1 
Faraday constant, F (C mol7!) 96,487 
GDL porosity 0.6 
Porosity of catalyst layer 0.12 
Volume fraction of ionomer in catalyst layer 0.4 
GDL permeability (m7) 1.0E—12 
Catalyst layer permeability (m°) 1.0E—13 
Equivalent weight of ionomer (kg mol!) 1.1 
Dry membrane density (kg m~?) 1980 
Effective electronic conductivity in CL/GDL (S m7!) 5000 
Operation pressure (atm) 2 
Condensation rate coefficient (s7!) [26] 5000 
Evaporation rate coefficient (s7! Pa~!) [26] 1.0E—4 
Liquid water density (kg m~’) 1000 
Liquid water viscosity (N s m~?) 3.5E—4 
Surface tension (N m~!) 6.25E—2 
Contact angle in GDL 110 
Contact angle in CL 95 
Thermal conductivity of GDL (W m7! K7!) [26] 1.5 
Thermal conductivity of CL (W m~! K!) [26] 1.5 
Thermal conductivity of the membrane (W m~! K7!) [26] 0.5 
Heat of vaporization (J kg~!) 2.3E+6 
Density of carbon material (kg m7) [32] 2200 
Heat capacity of carbon material (Ikg~! K7!) [32] 1050 
Heat capacity of liquid water (Ikg~! K7!) [32] 4200 
Heat capacity of membrane material (Jkg~! K7!) [32] 1050 

Electron transport: 

y eff ne 

‘(a Vds) + Ss = 0 (16) 


Based on the theoretical analyses in [27,33], since the electro- 
chemical double-layer charging/discharging process is very fast, 
the transient terms in the two equations can be safely neglected 
and the two equations thus remain in the steady-state form. 

The relevant expressions for the source terms and the other 
physicochemical relationships can be found in Table 2 and the 
reference [26]. 

The conservation equations, Eqs. (1)—(3), (11)-(13), (15) and 
(16), constitute the present transient two-phase non-isothermal 
PEM fuel cell model. This model is applied for numerical simu- 
lations in a two-dimensional configuration to fully examine the 
effects of liquid water transport and heat transfer phenomena on 
the transient responses of a PEM fuel cell. The computational 
domain includes five regions, namely the gas diffusion layers 
and catalyst layers on both the anode and cathode sides, and 
the membrane. In the present two-dimensional simulation, since 
fluid flows in gas channels are neglected, the gaseous velocity 
and pressure are not computed. This simplification is valid based 
on the theoretical and numerical analyses of the Peclet number 
in the porous materials in [13,37]. Therefore, the conservation 
equations, Eqs. (3), (11)—-(13), (15) and (16), are actually solved 
in the present two-dimensional simulations. 


Cathode 


Fig. 1. Cell geometry and the related boundaries. 


There are five boundary conditions to be specified, as shown 
in Fig. 1. The details regarding these boundary conditions can 
be found in [26] and are thus not repeated in this paper. 

The present numerical studies will focus on the transient 
responses of a PEM fuel cell with a step change of cell volt- 
age under various operation conditions, i.e. single-phase versus 
two-phase and isothermal versus non-isothermal. The initial 
conditions are defined as the steady-state numerical solutions 
at the previous cell voltage. 


3. Result and discussion 


The present transient two-phase non-isothermal PEM fuel 
cell model has been implemented into a commercial CFD pack- 
age, Fluent, through its user coding capabilities and applied 
herein for two-dimensional numerical simulations, as shown in 
Fig. 1. The geometric parameters of the fuel cell are listed in 
Table 3. 

Hydrogen and water vapor is fed into the anode while air and 
water vapor into the cathode. The fuel cell is operated at 2 atm 
on both the anode and cathode sides. The cell stoichiometry 
number is two on both sides with a reference current density of 
1 Acm~”. For all the calculations carried out in this paper, the 
boundary temperature is fixed at 80°C. In order to investigate 
the effects of two-phase flows and liquid water flooding on the 
transient responses of the fuel cell, the inlet gas mixtures are set 
as fully humidified at an inlet temperature of 80 °C on both sides. 
Based on these operation conditions, the inlet species concentra- 
tions can be easily determined and specified at boundaries | and 
2. Therefore, the present transient two-dimensional numerical 
calculations simulate a cross-section perpendicular to the flow 
direction at the inlet of a PEM fuel cell. 

Careful grid independence study has been conducted in [26], 
and a total of 1600 computational cells have been used in the 
present simulations. In order to obtain time-accurate solutions, 
variable time steps have been used in the present calculations. 
The time step is initially very small, 0.01 s, in order to cap- 
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Table 2 
Electrochemical and physical relationships 
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Description 


Expression 
H2 


ef 


Qa + Ae 


Transfer current density (A m7?) j= aloa 


1/2 
c] 2 
(=) ( 


ref 


RT 


F. n) , inanode side 


CO, ( 
- exp 
CO), ref 


ie F. n) ,  incathode side 


J = ajoe RT 
Over-potential (V) N = s — e, inanode side 
n = $s — e — Uo, in cathode side 
Open-circuit potential (V) Uo = 1.23 — 0.9 x 107°(T — 298) 
1.0, forà < 14 
Electro-osmotic drag coefficient na = = (Q —14)+1.0, otherwise 
CyRuT 
Water activity a= ahaa 
ps 
Water saturation pressure (atm) logio ps = —2.1794 + 0.02953(T — 273.15) — 9.1837 x 10->(T — 273.15)? + 1.4454 x 10-7(T — 273.15) 
Partial pressure of water vapor (Pa) pY = CwRaT 
3.1 x 1077A (928% — 1) - el -2346/T] 0<A<3 
Membrane water diffusivity (m? s~!) y= N 3 
4.17 x 1078A(1 + 161 e*) - el72346/T], otherwise 
Water content diffusivity (mol m`! s7!) D, = oy ee 


Proton conductivity (S m7!) 


k = (0.51394 — 0.326) exp [1268 ( 


1 L] 
303 T 


ture the initial rapid variations, and it increases gradually as the 
calculations progress. 

The transient responses of the PEM fuel cell corresponding to 
a step change of cell voltage are first studied under both single- 
phase and two-phase isothermal conditions. In the two-phase 
calculations, the liquid droplet coverage phenomenon at the gas 
diffusion layer (GDL) and gas channel (GC) interface is tem- 
porarily neglected, meaning that sint on the cathode side is set as 
zero. The effects of the liquid droplet coverage at the GDL/GC 
interface on the cathode side will be discussed in detail later in 
this section. 

Fig. 2 compares the dynamic responses of the fuel cell with 
a step decrease of cell voltage from 0.75 to 0.65 V under both 
single- and two-phase conditions. Since the contact resistances 
from the electron transport process are neglected in the calcula- 
tions, the cell performs much better than it could in real-world 
operations, but this should not affect the present numerical stud- 
ies. Results in Fig. 2 clearly indicate that under the two-phase 
condition, although the dynamic response time does not change, 


Table 3 
Cell geometric parameters 


Fuel cell geometry [mm] 


Layer thickness 
Diffusion 0.3 
Catalyst 0.01 
Membrane 0.025 
Land width 0.5 
Channel width 1.0 
Computational cell numbers ~1600 


the magnitude of the current density over-shoot increases, com- 
paring to that under the single-phase condition. As discussed 
in [27,32], the magnitude of the current density over-shoot 
results from the initial richer oxygen concentration maintained 
at a higher cell voltage of 0.75 V. Under the two-phase con- 
dition, with the decrease of the cell voltage and consequently 
the increase of the current density, the liquid saturation value 
increases as well, eventually rendering less oxygen reaching 
the catalyst layer at the steady-state and thus increasing current 
density over-shoot. 

Fig. 3 displays variations of the oxygen concentration in the 
cathode catalyst layer at different time steps under the two-phase 
condition. The results indicate that the oxygen concentration 
changes dramatically from 0 to 0.25s, but after that, it stabi- 


(a) F 
1.8 sH 
£ N 
oO 
2 17 F 
> L 
= L 
S H 
A 0.75V - 0.65V F 0.75V - 0.65V 
= 16 6F 
£ [ 
= 
S] b 
ES Sr 
-1 0 1 2 ae | 
Time(s) Time(s) 


Fig. 2. Transient responses of the fuel cell with a step decrease of cell voltage 
under (a) single-phase and (b) two-phase isothermal conditions. 
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0.25s 1.6s 


AQ 
Qa 
a 
w 


42 
X 
= 


Fig. 3. Transient variations of the oxygen concentration in the cathode catalyst 
layer corresponding to the same operation condition as in Fig. 2b. 


lizes inside the cathode catalyst layer. This clearly indicates that 
the oxygen concentration dictates the magnitude of the current 
density over-shoot but not the dynamic response time. In fact, 
the dynamic response time is controlled by the membrane hydra- 
tion/dehydration process, as manifested in Fig. 4, which displays 
the transient variations of the water content inside the mem- 
brane under the two-phase condition. With the decrease of the 
cell voltage from 0.75 to 0.65 V, the current density increases, 
and as such the electro-osmotic drag increases as well, causing 
larger water content gradient inside the membrane. The water 
content inside the membrane varies from 0 to 1.6s, consistent 
with the cell dynamic response time. 

Fig. 5 compares the dynamic responses of the fuel cell with 
a step increase of cell voltage from 0.65 to 0.75 V under both 
single- and two-phase conditions. Now the current density expe- 
riences under-shoot initially under both conditions, dictated by 
the low oxygen concentration initially maintained at the cell 
voltage of 0.65 V. Because of the presence of liquid water, the 
magnitude of the current density under-shoot is slightly stronger 
under the two-phase condition than that under the single-phase 
condition. In addition, the results in Figs. 2 and 5 indicate that 
the cell response time is less than 2s with a step change of 
cell voltage under both the single- and two-phase conditions. 
This further verifies that using a thin membrane would reduce 
the transient response time of a PEM fuel cell, a conclusion 
consistent with the results in reference [32] calculated using a 
single-phase model. 

The effect of heat transfer process on the transient response 
of a PEM fuel cell is further examined under the two-phase 


Os 0.25s 1.6s 


Qa 


Fig. 4. Transient variations of water content inside the membrane corresponding 
to the same operation condition as in Fig. 2b. 
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Fig. 5. Transient responses of the fuel cell with a step increase of cell voltage 
under (a) single-phase and (b) two-phase isothermal conditions. 


condition in Fig. 6, which compares the transient variations 
of the current density with a step decrease of cell voltage 
from 0.75 to 0.65 V under both isothermal and non-isothermal 
conditions. Although under the non-isothermal two-phase con- 
dition, the magnitude of the current density over-shoot only 
increases slightly, the response time of the fuel cell increases 
dramatically, from around 1.6s under the isothermal condition 
to around 8s under the non-isothermal condition. The same 
conclusion can be drawn from Fig. 7, which compares the tran- 
sient responses of the current density with a step increase of 
cell voltage from 0.65 to 0.75 V under both isothermal and 
non-isothermal two-phase conditions. In addition, under both 
isothermal and non-isothermal two-phase conditions, the tran- 
sient response time with a step decrease of cell voltage, as shown 
in Fig. 6, is longer than that with a step increase of cell voltage. 

It can be concluded from the results in Figs. 6 and 7 that 
heat transfer process and its intricate interactions with the liq- 
uid water saturation increase the transient response time of a 
PEM fuel cell undergoing a step change of cell voltage. This 
conclusion is consistent with that drawn in [32] calculated using 
a single-phase model. Fig. 8 illustrates the temperature distribu- 
tions at different time steps corresponding to the same operation 


(a) 1.9 (b) 1.9 


oe 
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Current Density(A/cnt ) 
G æ 
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Current Density(A/cm^) 
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a 
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0.75V - 0.65V 


0 2 4 6 8 10 


Time(s) Time(s) 


Fig. 6. Transient responses of the fuel cell with a step decrease of cell voltage 
under two-phase (a) isothermal and (b) non-isothermal conditions. 
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Fig. 7. Transient responses of the fuel cell with a step increase of cell voltage 
under two-phase (a) isothermal and (b) non-isothermal conditions. 


conditions as in Fig. 6b. Significant temperature variations can 
be observed from 0 to around 2.0 s on both the anode and cath- 
ode sides, consistent with the dramatic variations of the current 
density shown in Fig. 6b. After 2.0 s, only slight variations of 
the temperature distribution can be found, also consistent with 
the small current density variations in Fig. 6b. The temperature 
variations are closely coupled with the liquid water saturation 
variations, as shown in Fig. 9. With the temperature increasing, a 
condensation/evaporation interface appears in the porous mate- 
rials directly under the gas channel, and a dry region expands 
with time, especially from 0 to 2.0s. 

The transient variations of temperature and liquid satura- 
tion distributions would affect the water content distribution 
inside the membrane, as shown in Fig. 10. With the increase 
of the temperature and the appearance of a dry region directly 
under the gas channel, the water content inside the membrane 
decreases at the corresponding region, as clearly illustrated in 
Fig. 10. Therefore, it can be concluded that the heat trans- 
fer process influences the transient responses of a fuel cell 
through its intricate interactions with liquid water saturation 
in the porous materials and the water content inside the mem- 
brane under two-phase non-isothermal conditions, which in 
turn will affect the electrochemical kinetics, species trans- 
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357.2 
356.9 
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354.9 
354.6 
354.3 
354.0 


Fig. 8. Transient variations of cell temperature corresponding to the same oper- 
ation condition as in Fig. 6b. 


port, and membrane hydration/dehydration processes in a PEM 
fuel cell. 

The numerical results presented above are all calculated with- 
out considering the liquid droplet coverage phenomenon at the 
GDL/GC interface. In fact, downstream in the gas channel, with 
the water vapor oversaturated and the decrease of current density 
and consequently the decrease of temperature, the condensa- 
tion/evaporation interface directly under the gas channel would 
disappear, and as such liquid water will emerge from the GDL 
and attach at the GDL/GC interface. This phenomenon has been 
fully examined and discussed in reference [17] for a steady- 
state two-phase PEM fuel cell simulation. In order to include the 
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Fig. 9. Transient variations of liquid water saturation in cathode GDL corresponding to the same operation condition as in Fig. 6b. 
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Fig. 10. Transient variations of water content inside the membrane corresponding to the same operation condition as in Fig. 6b. 


effects of the liquid droplet coverage on the transient responses 
of a PEM fuel cell in the present two-dimensional simulations, 
the heat transfer process is neglected in the following calcula- 
tions. Fig. 11 compares the transient responses of a PEM fuel cell 
undergoing a step decrease of cell voltage from 0.75 to 0.65 V 
under the two-phase isothermal condition with and without the 
liquid droplet coverage phenomenon. The transient current den- 
sity variations with a step increase of cell voltage from 0.65 to 
0.75 V, with the other operation conditions remaining the same 
as in Fig. 11, are presented in Fig. 12. In both Figs. 11b and 12b, 
the liquid saturation values defined at the GDL/GC interface 
are 0.5 at 0.65 V and 0.3 at 0.75 V, appropriately accounting for 
the current density effect on the degree of the interfacial liquid 
droplet coverage. As the liquid droplet coverage at the GDL/GC 
interface strongly affects the liquid saturation value in the porous 
materials, it will significantly affect species transport, resulting 
in drastic over-shoot and under-shoot of the current density cor- 
responding to a step change of cell voltage. The over-shoot of the 
current density with a step decrease of cell voltage is especially 
strong, reaching around 20% of the final steady-state value, as 
shown in Fig. 11b. Both Figs. 11b and 12b indicate that the 
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Fig. 11. Transient responses of the fuel cell with a step decrease of cell voltage 
under isothermal two-phase conditions; (a) without interfacial liquid droplet 
coverage and (b) with interfacial liquid droplet coverage. 
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Fig. 12. Transient responses of the fuel cell with a step increase of cell voltage 
under isothermal two-phase conditions; (a) without interfacial liquid droplet 
coverage and (b) with interfacial liquid droplet coverage. 


effect of the interfacial liquid droplet coverage on the transient 
response time is negligible. However, it should be emphasized 
that in the present numerical simulations, the liquid saturation 
value at the GDL/GC interface has been assumed to change 
instantly with a step change of cell voltage. In the real-world 
operations, the interfacial saturation value should vary with time, 
and therefore, it could exert influences on not only the magnitude 
but also the time of the transient response process. The transient 
characteristics of the liquid droplet coverage at the GDL/GC 
interface should be determined through extensive experiments, 
and deserves further studies. Based on the present numerical 
results and discussions, it can be concluded that the effect of 
the liquid droplet coverage at the GDL/GC interface on the 
transient response of a PEM fuel cell is extremely strong and 
thus cannot be neglected in two-phase transient PEM fuel cell 
simulations. 


4. Conclusion 


In this paper, a transient two-phase non-isothermal PEM 
fuel cell model has been developed based on the previously 
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established two-phase mixed-domain approach. This model is 
capable of solving two-phase flow and heat transfer processes 
simultaneously with a proper consideration of the effect of liquid 
droplet coverage at the gas diffusion layer and the gas channel 
interface. The model has been applied for two-dimensional time- 
accurate simulations in a cross-section perpendicular to the flow 
direction to fully examine the effects of liquid water transport 
and heat transfer phenomena on the transient responses of a 
PEM fuel cell undergoing a step change of cell voltage, with 
and without condensation/evaporation interfaces. 

The present numerical results show that under isothermal 
two-phase conditions, the presence of liquid water in the porous 
materials increases the current density over-shoot and under- 
shoot, compared with the results under isothermal single-phase 
condition. Under the non-isothermal two-phase conditions, 
the heat transfer process significantly increases the transient 
response time. The heat transfer process influences the transient 
responses of a PEM fuel cell through its intricate interac- 
tions with liquid water saturation in the porous materials and 
the water content inside the membrane, which in turn affect 
the electrochemical kinetics, species transport, and membrane 
hydration/dehydration processes in the fuel cell. 

The present studies also indicate that, with an inclusion of the 
liquid droplet coverage at the GDL/GC interface, the interfacial 
liquid saturation level strongly affects the liquid saturation value 
inside the porous materials, and it in turn significantly influences 
species transport, resulting in drastic over-shoot and under-shoot 
of the current density with a step change of cell voltage. In fact, 
if the transient characteristics of the interfacial liquid droplet 
coverage could be determined and applied in the simulations, it 
would exert influences on not only the magnitude but also the 
time of the transient response process. This phenomenon will be 
further investigated in the future research work. 
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